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ABSTRACT 

Late asymptotic giant branch (AGB) and early post-AGB stars which are 
progenitors of planetary nebulae lose mass at extremely high rate, in what is 
termed a superwind. We show that the existence of this superwind during the 
post-AGB phase cannot be explained with models where the stellar effective tem- 
perature is the sole main physical parameter which determine the mass loss rate. 
Instead, we argue that the envelope structure, in particular the entropy and den- 
sity gradients, should be among the main parameters which determine the mass 
loss rate on the tip of the AGB and the early post-AGB evolutionary phases. 
The entropy profile becomes steeper and the density profile becomes shallower 
as the star becomes hotter on the early post-AGB phase, until the star heats-up 
to T > 8000 K. We propose that mass loss rate stays very high because of the 
envelope structure, and drops only when the effects of the temperature becomes 
important once again as the post-AGB star heats to ~ 6, 000 K. We do not 
propose a new mass loss mechanism, but rather mention several mechanisms by 
which these profiles may influence the mass loss rate within the popular mecha- 
nism for mass loss on the AGB, where pulsations coupled with radiation pressure 
on dust cause the high mass loss rate. 

Subject headings: stars: AGB and post-AGB - circumstellar matter - stars: 
mass loss - planetary nebulae: general 



1. INTRODUCTION 

Models for the acceleration of winds from asymptotic giant branch (AGB) stars which 
are based on radiation pressure on dust agree well with observations (e.g., Elitzur, Ivezic 
& Vinkovic 2002). The exact mass loss rate depends on the process which initiates the 
outflow (Elitzur et al. 2002). One of the most popular mechanism to initiate the outflow, 
which is based both on observations and theoretical studies, is stellar pulsation (e.g., Wood 
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1979; Bowen 1988 Hofner et al. 1996; Hofner 1999; Winters et al. 2000; Andersen, Hofner 
& Gautschy-Loidl 2002). In such models, it is found that the mass loss rate substantially 
increases as the star evolves along the AGB (e.g., Bedijn 1988; Bowen & Willson 1991; 
Hofner & Dorfi 1997; Wachter et al. 2002, hereafter WA02), with strong dependence on 
the stellar mass, effective temperature, and luminosity (or radius). This behavior along the 
AGB is in agreement with observations which require the mass loss rate to substantially 
increase near the termination of the AGB, in what was termed superwind (Renzini 1981). 
The strong dependence on the temperature (e.g., WA02) implies that as the star evolves 
along the post-AGB track and becomes hotter, at constant luminosity, the mass loss rate 
steeply decreases with time. However, observations of planetary nebulae (PNs) and stellar 
evolution calculations along the post-AGB require the high mass loss rate to continue during 
the post-AGB phase (e.g., Tylenda & Stasinska 1994), until the star heats- up to T ~ 5000 K 
(Blocker 1995). Indeed, the post-AGB evolution time obtained by WA02 is much too long 
(R. Szczerba 2002, private communication). We elaborate on this in Section 2. 

It seems that additional physical parameter(s) are needed to determine the mass loss rate 
(temperature still should be a factor as it influence dust formation). This is also supported by 
the observations that some elliptical PNs have spherical halo— which was presumably formed 
from the AGB 'regular' wind— while the inner dense region— which was probably formed 
from a superwind— is axisymmetrical. The change in mass loss geometry suggests that other 
factors are involved in determining the mass loss rate and geometry. Soker & Harpaz (1999, 
hereafter SH99) propose that the density and entropy profiles are the additional factors. 
Because of the high mass loss rate the envelope mass decreases very fast, such that even when 
the star starts its contraction on the post-AGB early phase, the average density of the stellar 
envelope decreases. This makes the density profile below the photosphere much shallower, 
and the entropy profile much steeper. SH99 suggest that one possibility by which these 
profiles enhance mass loss rate is that they enhance magnetic activity. Other possibilities 
are a switch in the pulsation mode, possibly to include non-radial pulsation modes, or a 
more violent convection, because of the steeper entropy profile. A more violent convective 
region may influence the pulsation and deformation of the surface of the AGB star, and 
cause shocks to propagate through its photosphere (Jacobs et al. 2000), possibly increasing 
the mass loss rate. These mechanisms require more investigation. 

SH99 concentrate on the onset of the superwind until the stellar radius shrinks to R ~ 
260i? Q and the temperature reaches T ~ 3300 K. In the present paper we show that the same 
arguments used by SH99 for the onset of the final intensive wind— the superwind— imply that 
the superwind ends when the effective temperature is around T ~ 5000 K. We emphasize 
that we do not propose a new mass loss mechanism. We accept that pulsations coupled with 
radiation pressure on dust is the mechanism for mass loss, and that stellar mass, effective 
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temperature, and radius play a significant role. However, we suggest (Sect. 3) that the 
envelope structure, in particular the entropy and density profiles, are important physical 
factors in determining the mass loss rate on the early post-AGB phase. 

2. PROBLEMS WITH STRONG TEMPERATURE DEPEND ANCE 

In the physical models cited above for mass loss from upper AGB stars (e.g., Bowen & 
Willson 1991; WA02), the mass loss rate depends mainly on the stellar effective temperature, 
T, the stellar luminosity, L (or stellar radius R) , and the total stellar mass M*. These 
models cannot reproduce the correct mass loss rate after the star evolves to the left on 
the HR diagram, i.e., when it becomes hotter at a constant luminosity. The reason is as 
follows. In order to get a very large increase in the mass loss rate as stars approach the tip 
of the AGB, as inferred from observations, the mass loss rate should strongly depend on the 
different variables, and in particular on the temperature. WA02 find the mass loss rate, in 
M Q yr _1 , to be 

log Mt = 8.86 - 1.95 log M*/M - piogT/K + 2.47 log L/L Q , (1) 

with P = 6.81. The subscript T stands for winds where the stellar effective temperature is 
the main parameter determining the mass loss rate. An AGB star starts its blueward motion 
(to the left on the HR diagram) when its envelope mass is M env ~ 0.1 — 0.4M Q (Schonberner 
1983; Soker 1992); stars with a more massive core turn blueward while retaining a more 
massive envelope. With the strong dependence on the temperature, the mass loss rate 
decreases very fast. The result is that the evolution time between the maximum mass loss 
rate epoch and the PN stage is much too long. 

We now show that the post-AGB evolution time according to the mass loss formula of 
WA02 is too long. In the model calculated by Soker (1992) the star reaches the tip of the 
AGB when its envelope mass is M env o — O.16M , then stays at that radius until the envelope 
mass decreases to M env i ~ 0.08M Q . After that point the envelope shrinks according to 

f M \°' 2 

R - m \0m^J jRq ' f ° r M co„>KT 5 M , (2) 

where M con is the mass in the convective region. To a good accuracy the equality M con ~ M env 
holds until the envelope mass reduces to M cnv ~ O.OO2M , the temperature rises to ~ 4000 K, 
and the stellar radius shrinks to R ~ 18O-R . The stellar luminosity in that model is 
L = 75OOL , hence the temperature is T = 27OO( J R/4OO J R )- 1/2 fT. During the later phase, 
the stellar mass does not change much, and stays in the range of 0.6 — O.68M . Substituting 
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the temperature, luminosity, and mass at the lower end in the expression for the mass loss 
rate of WA02 (eq. 1), and using equation 2, gives for the mass loss rate 



where (3 is defined in equation 1. The total time required to lose an envelope mass of M env i, 
during which approximation (2) holds, is given then by (assuming (3 < 10) 



Because (3 is close to 10, the factor 1/ (1—0.1/9) is large, and the evolution time is too long. For 
the parameters cited above with (3 = 6.81, it takes tti — 10 4 yrs to lose the last O.O8M . This 
evolution is too long, as the dynamical age of the dense shell of many PNs is much below that 
value (Sect. 1). Also, the maximum mass loss rate according to equation (1), when applied to 
the model from Soker (1992), is reached already when M env0 = 0.16M©. The evolution time 
from that point to the point when M envl = 0.08M© is r T0 = (0. 16- 0.08) /3 x 10~ 5 = 2700 yr. 
Overall, the PN stage will start ~ 12000 yrs after the superwind has started. This is a too 
long time when compared to most PNs. Accurate calculations by WA02 find this evolution 
time for the more massive cores, and a much longer evolution time from the AGB tip to the 
PN stage for lower mass stars. This is in contradiction with observations (R. Szczerba 2002, 
private communication). Note that Bowen & Willson (1991) get an ever increasing mass 
loss rate because they don't treat properly the AGB and post-AGB evolution, as they don't 
include the blueward stellar evolution part. 

Although crude, the treatment above has a robust conclusion: Models which attribute 
the final intensive wind on the AGB and and post-AGB (the superwind) mainly to the 
decreasing effective temperature of the star (or its increasing radius), are in contradiction 
with observations of PNs and proto-PNs, because these models predict a too long post-AGB 
evolution. 




(3) 
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3. THE DEPENDANCE OF MASS LOSS RATE ON THE ENVELOPE 

STRUCTURE 

3.1. The density profile 

Following SH99, we examine the ratio of the density at the photosphere, p p , to the 
average envelope density p a . 1 . The photospheric density is given by (Kippenhahn & Weigert 
1990) 

_ 2 pm H GM* 
Pp ~ 3~fc^^T' [ } 

where prrin is the mean mass per particle, ks is the Boltzmann constant, and k is the opacity. 
Substituting typical values for a post-AGB star in the last equation gives 

ft = 93 x ir '° ( io-'^g- ) ' (soSjk) 3 {wrf (a^e) g cm " (6) 

During the post-AGB phase, most of the envelope volume is convective. We are interested 
in the density in that region, hence we consider the mass of the convective region, M con . The 
average density in the convective region is given by 

* = i^s = 6 '° x 10 " (booTTk) 6 ( Toy fe) g cm " 3 (7) 

The ratio of photospheric to average density is then 

e> * o oi6 ( i 1 1 (-?—) 3 (^-) 1/2 (J±-) 

Pa V 10_3cm2 g _ v vsoook; \io 4 l q ) \o.6m q J \o.im ( 

At an effective temperature of T ~ 2800 K, appropriate for the tip of the AGB, the opacity 
is k ~ 3 x 10~ 4 . We find from equation (8) that at the tip of the AGB, when most of 
the envelope is convective and M env = M con , that p v j ' p a ~ O.3(M env /O.lM ) _1 . Of course, 
the average density must be larger than the photospheric density, hence the star starts to 
contracts, i.e. becomes hotter, before this stage is reached, i.e., when the envelope still has 
a significant amount of mass. 

During the early post AGB phase, when T < 4000 K, the envelope shrink according 
to equation (2) and the opacity (e.g., Alexander & Ferguson 1994) can be approximated by 



1 Note that the density scale in Figs. 1-5 of SH99 is too low by a factor of 10; the correct scale is displayed 
in their Fig. 6. 
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k ~ 4 x 10- 4 (T/3000fsT) 4 cm 2 g" 1 . Using these in equation (8) with L = 75OOL , gives for 
the density ratio (SH99) 

Pp ~ o.3 ( J^L-\ ~ 0.4 ( T ] for 3000 K < T < 4000 K. (9) 
p a \0.1MqJ V 3000 k / 

This ratio increases with mass loss, and implies that the density profile becomes shallower, 
almost flat, as indeed found in the numerical stellar models of SH99 for the early post-AGB 
phase. In SH99, the evolution of the envelope was followed at the end of the AGB phase, 
where the envelope mass decreased down to O.O15M (where the mass of the stellar core was 
kept constant at 0.6M Q ). 

In the present work we follow the evolution during the early post-AGB. The numerical 
methods used here are the same as those used by SH99, but the present model has somewhat 
higher luminosity, L = 9300L© instead of L — 8400L© in SH99 or L = 75OOL in Soker 
(1992). To demonstrate the same behavior as found for late AGB stars for these later stages, 
we plot in Figure 1 the density profiles at several post-AGB evolutionary points, where the 
envelope mass dropped down to the values of O.OOO2M . The envelope mass at the points 
displayed in the figure, has the values of 0.0015, 0.001, 0.0005, 0.0003, and O.OOO2M . 
Similar figure is plotted for the entropy profile in the same models in Figure 2 (see next 
section). 

An important point confronted in such evolutionary calculations, concerns the fact that 
due to the very low density in the outer envelope, the opacity in this region is very low, 
and the effective temperature of the star, calculated at the point where the optical depth, 
r, reaches the value of 2/3, is deep in the outer envelope. In the figures, we display the run 
of the density and of the entropy out to the point where they actually do not play any role 
in the construction of the envelope structure. This point is about 5 — 15R Q further than the 
point where the effective temperature was calculated. 

An interesting feature that appears in the density run of all the models is the density 
inversion, that takes place in the outer part of the envelope. This phenomenon was discussed 
by Harpaz (1984). Close to the outer region where convection takes place, where the density 
is very low, the convective flux becomes very uneflicient, because it depends on the density 
(see eq. 10 below). A very steep temperature gradient develops in order to increase the 
convective flux. However, this steep temperature gradient, causes an increase of the pressure, 
and since the pressure gradient have to satisfy the hydrostatic equation, a density inversion 
is developed, in order to compensate for the increase of the pressure due to the increase of the 
temperature. Certainly, the demand on the increased temperature gradient is removed when 
the opacity falls steeply at the point where the recombination of the hydrogen is completed, 
and there, the density inversion vanishes. 
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At the temperature range of 5000 K < T < 6000 K, a steep increase of the opacity 
implies that the photospheric density steeply decreases with increasing effective temperature. 
Using the opacity table of Alexander & Ferguson (1994), we find from equation (6), with the 
same luminosity and stellar mass used there, that at T = 4000 K, k = 10~ 3 cm 2 g _1 and 
p p ~ 9 x 10~ 10 g cm~ 3 , at T = 5000 K, « = 2.2 x 10~ 3 cm 2 g" 1 and p p ~ 4 x 10~ 10 g cm" 3 , 
and at T = 6000 K, k = 0.03 cm 2 g _1 and p p = 3 x 10~ n g cm -3 . This demonstrates 
the sharp drop in the photospheric density occurring as the star moves blueward across 
T ~ 5000 K. Comparing the change in photospheric density with the average density in the 
convective region (Fig. 1), shows that the photospheric density drops faster than the average 
density in the convective region, as the post-AGB star heats up in the temperature range 
of T ~ 5000 — 6000 K. Over all, in evolving from the AGB tip blueward to a temperature 
of T ~ 5000 K the general envelope density profile becomes shallower, as evident from the 
increase of the ratio p v j p a . In the evolution from T ~ 5000 K to Tf ~ 7000 — 9000 K, with 
higher temperatures for more luminous post-AGB stars, the trend of this behavior changes, 
but the density profile in the convective region stays shallow (Fig. 1). Becoming hotter than 
Tf, the star contracts very fast and the envelope density becomes very steep (Soker 1992). 
The range relevant to the present paper is the post-AGB track cooler than Tf. 

If the shallow density profile is behind the high mass loss rate during the late AGB and 
early post-AGB phase (SH99; sec. 1 above), then it is expected that the mass loss rate will 
decrease significantly only as the star gets hotter than T = Tf ~ 7000 — 9000 K. Including 
the effect of temperature-low temperature is required to form dust-the sharp drop in the 
mass loss rate can occur somewhat earlier, at T ~ 6000 K. 



3.2. The entropy profile 

SH99 show that during the late-AGB and early post-AGB phases the entropy gradient 
becomes steeper with mass loss. In Figure 2 here we show that the entropy profile stays 
steep for most of the relevant post-AGB phase. 

For stars on the upper AGB and early post-AGB phases the convective velocity, v c , in 
the mixing length prescription is limited to not exceed the isothermal sound speed, and the 
pressure profile is relatively shallow (e.g., figs. 1-5 in SH99). Under these circumstances, the 
convective luminosity at each radius r in the envelope can be written as 

£co„v(r) ^ 4irar 2 T e (r)v c (r)p e (r)H P (r)\dS(r)/dr\, (10) 

where T e and p e are the temperature and density, respectively, inside the envelope (rather 
than on the photosphere), H P is the pressure scale height, and a is a parameter of order 



- 8- 



unity which is basically the ratio of the mixing length to Hp. The entropy S is in units 
of erg K" 1 g -1 . The envelope temperature profile, T e (r), during the relevant post-AGB 
evolution segment and in the relevant envelope parts, evolves very slowly (e.g., Soker 1992). 
The same holds for Hp and v c (SH99). Since the density in the envelope is very low, equation 
(10) implies that the entropy profile must be steep, as indeed we find from accurate numerical 
calculations (Fig. 2). 

Again, If the steep entropy gradient is an important factor in creating the high mass 
loss rate during the late AGB and early post-AGB phase (SH99; sec. 1 above), then it is 
expected that the mass loss rate will decrease significantly only as the star gets hotter than 
T = Tf ~ 7000 — 9000 K (higher temperatures for more luminous stars). As mentioned in 
the previous subsection, effects of dust formation due to effective temperature will cause a 
sharp drop in the mass loss rate already at somewhat lower temperatures. 

4. SUMMARY 

The main goal of the present paper is to argue that the decreasing temperature and 
increasing radius can't be the only main parameters which determine the final intensive 
wind— the superwind— on the tip of the AGB and the early post-AGB phases. The decreasing 
temperature and increasing radius, as well as the decrease in total stellar mass, lead to an 
increase in the mass loss rate of red giant branch and AGB stars. We accept this general 
trends. However, we showed (Sect. 2) that models with mass loss rates which depend mainly 
on these variables predict that the mass loss rate steeply decreases as the star leaves the AGB 
tip and becomes hotter (at constant luminosity). Because at that stage the envelope still 
contains a significant mass, the post-AGB evolution predicted by such models (e.g., WA02) 
is much too long. From observations and stellar evolution models fitted to observations, it is 
known that the final intensive wind should continue well into the post-AGB phase (Tylenda 
& Stasinska 1994), until the stellar effective temperature becomes T ~ 5000 — 6000 K, (higher 
temperature for more massive stars, e.g., Blocker 1995). 

Following SH99, we showed in this paper (Sect. 3) that the density profile stays shal- 
low, and the entropy profile stays steep until the post-AGB stars heats up above T = Tf ~ 
7000 — 9000, when the star shrinks fast (Soker 1992). Therefore, we suggest that the enve- 
lope structure, in particular the entropy and density gradients, should be among the main 
parameters which determines the mass loss rate on the tip of the AGB and the early post 
AGB. Some mechanisms which may influence the mass loss rate, and which are influenced 
by the envelope structure were mentioned in Section 1. 
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The sharp drop in the mass loss rate is likely to occur earlier, before the star heats-up 
to T = Tf, since it becomes more difficult to form dust as the star gets hotter. Over all, 
we suggest a general form for the mass loss rate on the tip of the AGB and early post-AGB 
phases 

M = M T [l + G(S',p% (11) 

where is the mass loss rate mainly determined from the effects of the effective temper- 
ature, stellar radius and mass (e.g., eq. 1 above taken from WA02). G(S',p') is a function 
of the density and entropy gradients in the convective envelope, with G ^> 1 only when the 
density profile is shallow and the entropy profile is steep in the envelope convective region. 
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Fig. I. — The run of the density in the envelope along five evolutionary points of a post-AGB star. 
The core mass in all the models is O.6M , the stellar luminosity is L = 93OOL , and the envelope 
mass at each evolutionary point is marked on the graphs. The density is plotted at several solar 
radii above the photosphere. Note the very sallow density profile in a large fraction of the envelope, 
which, we argue, influences the mass loss process. 
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Fig. 2. — Like Figure 1 but for the entropy instead of the density. Note the relative steep entropy 
profile, which, as argued in the text, influences the mass loss process. 



